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Charmonium dissociation cross sections due to flavor-exchange charmonium-baryon scattering 
are computed in the constituent quark model. We present results for inelastic J/^N and r) c N 
scattering amplitudes and cross sections into 46 final channels, including final states composed of 
various combinations of D, D* , E c , and A c . These results are relevant to experimental searches for 
the deconfined phase of quark matter, and may be useful in identifying the contribution of initial cc 
production to the open-charm final states observed at RHIC through the characteristic flavor ratios 
of certain channels. These results are also of interest to possible charmonium-nucleon bound states. 
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I. INTRODUCTION 



The production of heavy quarkonium in heavy ion collisions has long been considered a possible diagnostic for the 
appearance of exotic QCD phases in relativistic heavy ion collisions [HQ- I n particular it has been anticipated that 
' the charmonium production cross section in central AA reactions will be suppressed if a quark gluon plasma (QGP) is 
formed, since the long-range cc confining potential will be screened within the QGP. However, this diagnostic can be 
t-H ■ confounded by subsequent charmonium dissociation due to inelastic hadron rescattering "comover absorption" , which 
also contributes to the depletion of initially produced charmonium. Alternatively, charmonium can be regenerated 
due to rescattering of open-charm hadrons in the late stages of heavy ion collisions [3|. It is clear that a thorough 
understanding of these soft hadronic final state interactions is required before one can make confident statements 
regarding QGP production based on charmonium production cross sections [J]. The usefulness of an understanding 
of these soft processes is now expected to extend well into the deconfined phase, since recent lattice computations 
, indicate that low lying charmonia survive as resonances up to T « 3T C Q ■ 

Unfortunately, little experimental information exists regarding these charmonium dissociation cross sections. Some 
simple phenomenological estimates based on 'pre- vector meson dominance' @ and absorption cross sections in heavy 
ion collisions Q give afj^ oc s 22 mb and af^{y/s « 10 GeV) « 6 mb respectively. Alternatively, J/ip photoproduc- 
tion yields a J /tp N cross section of approximately 3.5 mbQ while a combined analysis of J/ip production from p + A 
collisions gives a result of 7 mb @. For more detailed predictions of the relevant near-threshold cross sections one 
must employ theoretical models of these scattering processes. 

Theoretical estimates have employed a variety of methods, and (perhaps not surprisingly in view of the lack of low- 
energy experimental data) predict cross sections that vary over several orders of magnitude in the relevant kinematic 
regime. Early estimates by Kharzeev and Satz [11] using the color-dipole diffractive model of Bhanot and Peskin 
|lfj | (which is only justified at high energies) gave extremely small near-threshold J/ipN total cross sections, typically 
of microbarn scale. More recently, JjipN dissociation cross sections have been estimated using meson exchange 
models [l2l . [l3j], assuming for example ^-channel charmed meson exchange or an SU(4)-symmetric hadron effective 
lagrangian. Results from these models for meson- J scattering cross sections near thresholds are typically in the few 
mb range. Similar computations have been reported for J/ipN scattering [l3l fl4T|. which also find total cross sections 
near threshold in the few mb scale. Although these meson exchange models are of great interest as possibly realistic 
descriptions of these near-threshold processes they suffer from uncertainties due to poorly understood vertex form 
factors, Fock space truncations in the set of exchanged particles, and the questionable assumption of higher symmetry 
groups such as SU(4) (this implicitly assumes a close relation ship between the dynamics of pions and heavy-quark 
cc mesons). Regarding hadronic form factors, several groups |16| have found a strong suppression of the predicted 
dissociation cross sections on incorporating plausible hadronic form factors. 

Many of the problems encountered in previous approaches are avoided if one implements a "microscopic" quark- 
gluon description of these scattering processes, for example using the constituent quark model. The earliest application 
of this approach to charmonium dissociation is the work of Martin et aLpjJ, who applied the method of Ref.[H[ to 
J/^7r — > D*D and D*D* scattering. Martin et al. assumed that the confining interaction only operated between 
q-q pairs; a more conventional quark-gluon model of charmonium dissociation cross sections based on the usual A • A 
color structure has been developed by Wong et al. These approaches describe hadronic interactions and bound 

states in terms of the nonrelativistic quark model, and usually assume that the scattering amplitudes are given to 
sufficient accuracy at Born order in the quark-gluon interaction (this can be relaxed of course). At Born order 
in these valence qq annihilation free channels hadron-hadron scattering occurs by constituent interchange, and the 
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scattering amplitudes and cross sections can be derived analytically if sufficiently simple wavefunctions are employed 
(for example simple harmonic oscillator forms). Applications of this approach to a wide range of similar scattering 
processes without valence annihilation, such as 1=2 irn [H, 22], KN [2(| and NN [2l[ have shown that this leads to 
numerically realistic results for many short-ranged S-wave scattering processes. 

In this paper we report the first results for charmonium dissociation due to scattering from nucleons in this type of 
constituent quark scattering model [23j . These amplitudes are also relevant to the time-reversed process of charmonium 
regeneration, and may prove useful in future studies of bound states of charm and nuclei, as suggested by Brodsky et 
al. [p. 

II. QUARK MODEL SCATTERING FORMALISM 

The Born-order quark interchange model approximates hadron-hadron scattering as due to a single interaction of 
the standard quark-model interaction Hamiltonian Hj between all constituent pairs in different hadrons [l8j . For the 
calculations reported here we employ a constituent quark model interaction of the form 

= E (~ - %i - o^f-e-^ Sl • Sj ) Ti • Tj (1) 
y V '•..» 4 :i\ 77111,111, J J J 

where the sum extends over all quarks and antiquarks. The three terms in this expression are respectively the 
color Coulomb interaction, linear confinement, and a regularized contact spin-spin hyperfine interaction. The model 
parameters have elsewhere been fitted to meson and baryon spectra (see Appendix lAl for their values). 

The color structure of Eq. [T]is given by the usual quadratic perturbative form T ■ T, where T = A/2 for a quark 
and — A*/2 for an antiquark (A is a Gell-Mann matrix). Meson-baryon systems do not have trivial color dynamics 
(i.e., a fixed color state with an overall constant color factor), in contrast to mesons and baryons individually, so 
the assumed color structure of EqQ] affects the relative amplitudes for different channels. Of course this T ■ T color 
structure represents a severe truncation of the full dynamics of the gluonic degrees of freedom. This model does 
however reproduce the relevant low-energy features of more complete models, and the T ■ T form is known to be 
realistic in describing lattice confinement potentials as well as the interactions of low lying hadrons. 

We now consider the generic inelastic two-body scattering process AB — ► CD, where A is a cc charmonium state, 
B is a nucleon, C is an open-charm (nc) meson, and D is an open-charm (nnc) baryon. (Here, n is a light u or 
d quark.) It is convenient to label the quarks as [cc] 12 [qi qi 93] 345. We choose the (45) quark pair to have definite 
symmetry under P45 quark interchange. Quark interchange scattering in all the systems considered here involves the 
specific quark permutations P23 and P23P34 = ^34^24- By inspection, P23 permutations give rise to P - (*'£+ + final 
states in J/ip(r) c )p collisions, and P/°^*- ) A°(E°) states in J/ip(r) c )n collisions. Similarly, J/ip(rj c ) interactions involving 
the permutation operator P23P34 give P°^*- ) A+(E+) final states, and from J/tfi(r] c )n one produces P~(*)A+(£+) final 
states. 

Due to P45 symmetry only four unique spatial matrix elements are encountered; these are summarized in Table HI 
The first row gives labels specifying the interaction potentials Vij . Thus for example d\ represents a spatial integral 
of the form J ipl 3 ip% 45 143^12^345 • 

TABLE I: Spatial Symmetries 



perm 


13 


14 


15 


23 


24 


25 


P23 




d 2 


<h 


di 


dz 


d 2 


P23P34 


<k 


di 


<h 


d-z 


di 


d 2 



The four unique matrix elements correspond to the four diagrams shown in Fig[TJ Note that quark line rearrange- 
ment is required to give a nonzero Born-order scattering amplitude, due to the color structure assumed for the quark 
model interaction. We follow the procedures described in Ref.[l8[ in evaluating these diagrams. We assume the 'prior' 
form of the scattering amplitude here, in which the separation of the full Hamiltonian into free and interaction parts 
is specified by the initial hadrons. 

Hadronic wavefunctions may gcncrically be written as linear combinations of product basis vectors of the form 
^ = CxS^, where the individual factors are the color, spin, flavor, and spatial wavefunctions respectively. Explicit 
wavefunctions for the mesons and baryons considered here are given in Appendix [5] Because of this factorizability 
of the hadron wavefunctions, which involve single factored terms in the cases we consider, the contribution of each 
quark diagram to the full scattering amplitude can be written as the sum of products of individual color, spin, flavor, 
and spatial matrix elements. 
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d\ d 2 
FIG. 1: The four quark-interchange meson-baryon scattering diagrams in the prior formalism. 

The Born order scattering amplitudes are then constructed by summing matrix elements of the interaction potential. 
For the case of P23 scattering, for example, we have 



A = -(P 23 E C E D \E A E B ) {P23C c C D \T l -T J \C A C B 

{P23XCXD\S 1 ■ Sj\XAXB) ■ (P2Slpci>D\Vij\ll> A ll>B 



(2) 



The spin factor Si ■ Sj is replaced by unity when considering central (pure potential) interactions. This expression 

can be simplified using the symmetries and notation discussed above to obtain A = —w ■ d, where d = (di, d 2 , d\, d 2 ) 
and w is a weight vector arising from flavor, color, and spin matrix elements. The weight vectors used in this work 
are listed in Appendix [B] 

Finally, the charmonium dissociation cross sections were computed from these amplitudes using the expression 



47T 2 k 



(3) 



where ^ab = E a Eb/{E a + Eb)- Charmonium regeneration cross sections can be also obtained using this formalism, 
since they are related to the dissociation processes by time reversal. 

III. CHARMONIUM DISSOCIATION CROSS SECTIONS 

A total of 46 exclusive charmonium-nucleon inelastic scattering processes were considered, involving nucleons and 
the J/tp or rj c in the initial state and all kinematically accessible open-charm S-wave-meson-baryon channels in the 
final state. Specifically, we have computed total cross sections for the following reactions: 



J/#; r)cP D°A+; L>°£+; L>°*A+; L>°*£+; D~*^ 



(4) 



J/tpn; r\ c n 



D°£°; £>°*£°; D"£+; £>"*£+; D~A+: L>~*A+ 



(5) 



1-V++ 



3/2 -> D *S c3/2 ;L> *S c3/2 ;D S c3/2 ;£> £ c3/2 



(6) 
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[J/ipr 



J3/2 



U lj c3/2> 1J ^3/2' U 3/2' U ^3/2 



(7) 



^c3/2 



and 



3/2- 



(8) 



(9) 



Isospin symmetry and the equivalence of many quark line diagrams imply that all the J/ip(ri c )n amplitudes are 
simply related to J/ip(r) c )p amplitudes, as follows: 



and 



A{J/^{r, c )n -> D°W£°) = A{J/^{r, c )p -> D"We++) 



A(J/iKVc)n - D-WE+(A+)) = A{J/^{r, e )p - S°«S+(A+)). 



(10) 



(11) 



Thus 23 unique amplitudes remain to be computed. 

An additional isospin relation between these reactions can be derived because the D~Y<+ + and -D°E^ states both 
couple to \II Z 



III)- Th us 



(12) 



This relationship evidently holds for the weights reported in Appendix [B] which were derived without making this 
assumption. 

The cross sections we find using this approach typically rise rapidly at threshold, and then are strongly damped 
above a hadron momentum scale set by the hadronic wavefunctions. Since the hypcrfinc, Coulomb, and linear 
confinement contributions vary in sign within channels and have very different momentum dependences, secondary 
peaks can appear (see Figs l2l3l) . although these tend to be much weaker than the near-threshold maxima. In three 
instances this pattern is reversed; one of these is shown in the right panel of FigfSJ 





E (Cev) 



E (Gev) 



FIG. 2: J/ipp -> D°A+ (left) J/ipp -> D°*A+ (right). Curves are: total cross section (solid), hyperfine (dotted), linear (dashed), 
Coulomb (small dash). 



We summarize all 23 independent cross sections in Tables IrflllVI The columns in these tables specify the energies 
and values of the near-threshold maxima (point 1), the subsequent minimum (point 2), and the secondary maximum 
(point 3), if it is significant. 
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FIG. 3: J/ipp -> D°*E+ (left) »? c p -> D _ E+ + (right). Curves 
Coulomb (small dash). 
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Final State 



D u At 
D°E+ 

D Q *A+ 
D°*E+ 
D°*E 
D 
D 



TABLE II: 
{E{GeV),a{mb))i 



r/ c p Cross Sections. 
{E(GeV),a{mb)) 2 (E(GeV),a{mb)) 3 



c3/2 

E++ 



e++ 



(4.196 , 
(4.355 , 
(4.360 , 
(4.311 , 
(4.5021 
(4.547 , 
(4.5050 
(4.551 , 



1.190) 
0.221) 
0.432) 
0.233) 
, 0.0073) 
0.018) 
, 0.0144) 
0.009) 



(4.527 , 0.006) 
(4.8162 , 0.0001) 
(4.702 , 0.003) 
(4.8195 , 0.0001) 
(4.706 , 0.001) 



(4.820 , 0.025) 
(5.1164 , 0.0003) 
(4.966 , 0.008) 
(5.1199 , 0.0005) 
(4.971 , 0.004) 



IV. DISCUSSION AND CONCLUSIONS 

The most important process is J/ipp — ► -D°A+ scattering, with a peak cross section of approximately 7 mb (Fig. 
[2]). Next in strength is J/ipp — > Z?°*A j~, w ith a maximum near 1.7 mb. These results disagree with meson exchange 
models, for example Liu, Ko, and Lin[14j find a larger cross section for D*A C scattering than DA C while Sibirtsev, 
Tsushima, and Thomas 14] find roughly comparable cross sections. 

The total dissociation cross section in shown in Fig. 2] It is clear that the quark model description of charm 
dissociation is not in agreement with the widely held view that the absorption cross section should be constant. In 
general cross sections increase with phase space and then are strongly damped due to the finite size of hadronic 
wavefunctions. Although this effect can be overcome by the opening of many channels, it is evident that this does not 
occur up to center of mass energies of 5 GeV or more. This observation is in disagreement with typical meson exchange 
models, which tend to find slowly increasing total cross sections. Presumably this behaviour is due to (approximately) 
cancelling momentum dependence between vertices with derivative coupling and monopole form factors. While the 
momentum dependence of quark model vertices must agree with their hadronic analogues (see, for example, Ref. [Hj]), 
quark model form factors fall much more rapidly in momentum than the power law monopole form factors, providing 
a possible source of the discrepancy. 

We remark that the quark exchange computation presented here can be regarded as the short range component of 
J /if) dissociation in scattering from nucleons. Long range dissociation can occur via virtual pion emission from the 
nucleon giving rise to reactions such as J/ipN — > DD*N . This amplitude can be related to that of J/ipn — > DD* 
scattering, which is relevant to comover absorption. This has been done, for example, by Liu, Ko, and Lin[l4j who 
obtain asymptotic cross sections of order 1 mb for each of J/tpN — ► DD*N, D*DN, and DDN. Of course, their results 
depend on theoretical J /ipit scattering cross sections as computed in the meson exchange formalism. Unfortunately, 
these are strongly dependent on form factors; for example J/ipn — > DD* ranges from 20 mb at large y/s with no form 
factor to 4 mb with a form factor cutoff of 1 GeV. Alternatively, quark model computations of the same reaction 
give maximum cross sections of order 1/2 mb just above threshold [19], in keeping with phenomenological estimates of 
comover absorption 0, [26[. From these considerations we conclude that the long range contribution to charmonium 
dissociation by nucleons is small (less that 1 mb). 

Comparison with experiment is difficult at present. As mentioned in the Introduction, J /ip photoproduction and 
p + A collisions yield estimates of 3.5 mb and 7 mb respectively. However these estimates are made assuming constant 
cross sections, and it appears unlikely that this is accurate. Nevertheless, the average near-threshold strength of our 
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Final State 



D°*At 
S°*E+ 
n°*Y+ 

U lj c3/2 

d-*e++ 



TABLE III: J/tpp (S = 1/2) Cross Sections. 
{E{GeV),(j{mb))i (E(GeV), a(mb)) 2 {E{GeV),a{mb))z 



(4.184 
(4.350 
(4.355 
(4.308 
(4.498 
(4.545 
(4.501 
(4.549 



6.815) 
0.126) 
0.246) 
0.120) 
0.180) 
0.009) 
0.354) 
0.017) 



(4.619 , 0.010) (4.884 , 0.022) 



(4.425 , 0.001) 
(4.668 , 0.002) 
(4.672 , 0.004) 



(4.703 , 0.079) 
(4.891 , 0.005) 
(4.896 , 0.009) 



Final State 



TABLE IV: J/ipp (S = 3/2) Cross Sections. 
(E(GeV),a(mb))i (E(GeV),a(mb)) 2 (E(GeV),a(mb)) 3 



3/2 
'-+ 
3/2 



D°*Af 
D°*E+ 

n°*y + 

U L c3/2 

n-*Y++ 

U ^3/2 



(4.418 
(4.424 
(4.324 
(4.498 
(4.545 
(4.501 
(4.549 



0.276) 
0.536) 
1.712) 
0.054) 
0.042) 
0.106) 
0.082) 



(4.873 , 0.000) 
(4.881 , 0.001) 
(4.639 , 0.021) 

(4.733 , 0.002) 
(4.737 , 0.004) 



(5.165 , 0.001) 
(5.173 , 0.002) 
(4.791 , 0.027) 

(4.986 , 0.006) 
(4.991 , 0.012) 



total cross section seems in accord with these phenomenological estimates. 

Finally, the strong near-threshold behaviour seen here implies that it would be worthwhile to revisit hydrodynamic 
simulations of charmonium suppression in heavy ion collisions. In this regard, we note that the current computation 
deals with free charmonium scattering, rather than the "pre" -charmonium expected in RHIC collisions. Further- 
more, substantial density dependence may be present [27| and this should be taken into account when applying these 
predictions to a nuclear environment. 
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APPENDIX A: HADRONIC WAVEFUNCTIONS 

Total baryon wavefunctions were constructed as follows: 



N = CZ N <S> BX \ (Al) 

A c = CZ A <5> DX P , (A2) 

£ c =CS E $r>x\ (A3) 

S c3/2 =CS E $ DX 3/2 (A4) 



where 
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Sa = —=c{ud — du), (A5) 
V2 



^ s ++ = cmt, (A6) 
S s + = -^c(ud + <iw), (A7) 
S s o = cdd. (A8) 



The baryon color wavefunction is given by 



C = -±=e abc (A9) 
V6 



where e is the rank three antisymmetric Cartesian tensor in color indices. 
Baryon spin wavefunctions are (only top states are given) 



X A n = + +-> + ! + -+> - 2| - ++)), (A10) 



2 2 



Xii = 4(1 ++->-! +-+», (AH) 



2 2 



and 



V5 V 

x| 7 a =! + ++>• (A12) 



The spatial wavefunctions are 



•*-#v*-HwtJ (A13> 

Q3/4 / „2 2 \ 

$B = 3/2 3/2 3/2 ^ -gf - S ( A16 ) 



with 



Prei = ; , (A17) 

mi + m,2 

_ /T (m 3 + 2m 5 ) p 4 - (m 3 + 2m 4 ) Ps + (ra 5 - m 4 ) p 3 fA1Sl 

Pp ~~ V 2 m 4 + m 5 + m 3 ' ^ J 

[3m 3 (p4+p 5 ) - (m 4 + m 5 )p 3 

P\ = \h; ; ; • A1 9 

V 2 m 4 + mg + m 3 

SHO meson and baryon wavefunctions can be expected to be reasonably accurate for hadrons composed of light 
quarks or heavy-light mesons. However, they are not reliable for Coulombic systems such as charmonium. Thus in our 
computations the nominal charmonium wavefunction ($^4 above) has been replaced with a sum over Gaussians that 
accurately replicates the numerically obtained wavefunction. The subsequent change in the cross section is illustrated 
in Fig. [5l where the importance of an accurate charmonium wavefunction near scattering threshold is evident. 

All model parameters have been fixed by meson and baryon spectroscopy. The parameters used here are listed in 
Table El 
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FIG. 5: a(ipp — * D Ac). Single J/ip SHO (dashed), six Gaussians (dotted), ten Gaussians (solid). 

TABLE V: Model Parameters. All units are GeV except b (GeV 2 ) and a s . 
m u m c b a a ffg I3 V [3 D fi D * ct N g p g x 
0.33 1.5 0.162 0.594 0.67 0.74 0.44 0.38 0.3 0.3 0.45 

Note that in an SHO (Isgur-Karl) model one has the relationships 

OL p = (3km) 1 / 4 

and 

a x = (3fcm A ) 1/4 

with 

3mM 

mx = 

2m + M 

where M is the mass of quark 3. 

APPENDIX B: WEIGHTS 

The color matrix elements required in the evaluation of the scattering amplitudes we evaluate arc listed in Table IVT1 
below. 

TABLE VI: Color Matrix Elements 
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With the spatial symmetries specified in Sect. [TT1 weights are computed as follows 

Wi = (P 2 3ScSd|SaS b ) (P23XcXd\xaxb) {P2zCcCd\T 2 ■ T 3 \CaCb) (Bl) 

and 

w 2 = (P 2 3ScSd|5a5 b ) (P23XcXd\xaXb) (P23CcCd\T 2 ■ Tk + T 2 • T b \C A C B )- (B2) 

Similar expressions hold for W3 and W4,. Here we consider the case of P23 scattering through a central potential. 
Substituting the color matrix elements of Table IVII then yields 



10 



w = -(P 23 E C E D \E A E B ) (P 23 xcXd\xaxb) (1, -1, -1, 1). (B3) 

The prefactor in this expression is reported in the column labelled 'Coulomb and linear' in the tables below. The spin- 
dependence of the hyperfine interaction does not permit simplification, hence the entire weight vectors are specified. 



TABLE VII: J/tpp, S = 1/2 Weights. 
Final State hyperfine Coulomb and linear 





$[-1,0,3,0] 


s3 


108 


[1,-2,-3,-2] 


h 
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% 
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V2 
108 


[7, 10, -5, -6] 


5V2 
27 
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2V3 

2 fo 


1 

54 


[-7,-10,5,6] 




i [1,1, 1,3] 


27 




-#[1,1,1,3] 


4s/2 
27 



D°A+ 
D°E+ 
D°*At 
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U ^c3/2 

r>-*y++ 



TABLE VIII: J/tpp, S = 3/2 Weights. 
Final State hyperfine Coulomb and linear 



"¥[1,0,-1,0] 
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-i [5,2,-1,0] 
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27 
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U ^3/2 
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V6 
54 
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Final State 



TABLE IX: rj c p Weights. 

hyperfine Coulomb and linear 
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